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Abstract 

One particle inclusive cross sections in the target fragmentation region are con- 
sidered and an explicit calculation is performed in ((/'^)6 model field theory. The 
collinear divergences can be correctly absorbed into a parton density and a fragmen- 
tation function but the renormalized cross section gets a large logarithmic correction 
as expected in a two scale regime. We find that the coefficient of such a correction 
is precisely the scalar DGLAP kernel. Furthermore the consistency of this result 
with an extended factorization hypothesis is investigated. 
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1 Introduction 



Asymptotic freedom has led to important predictions for hard inclusive hadronic pro- 
cesses in perturbative QCD. The basic tools are factorization theorems together with 
perturbative evaluation of scaling properties. Nevertheless perturbative QCD has been 
successfully applied also to various semi- inclusive processes, for which a complete factor- 
ization proof is still lacking. 

In a deep inelastic semi-inclusive reaction a hadron with momentum p is scattered by 
a far off shell photon with momentum q {Q^ = —q^ ^ ^qcd) ^^"^ from the inclusive final 
state, a hadron with momentum p' is detected. Until some years ago the main interest 
has been devoted in hadron production in the current fragmentation region, that is in the 
region in which the momentum transfer t = —{p — p'Y is of order 0. In the last few 
years, after the introduction of fracture functions 0, a novel attention has been payed 
to the target fragmentation region {t ^ Q'^) ^, With the observation of diffractive 
deep inelastic events at HERA part of this attention has been focused on the limit in 
which the longitudinal momentum fraction of the observed hadron is close to unity and t 
is of order 1 GeV^ or less. 

In this paper we want to deal with semi-inclusive deep inelastic scattering (DIS) in 
a rather simplified context, that is we will work in (0^)6- The scalar (0^)6 model share 
with QCD some important aspect. It is an asymptotically free theory and the Feynman 
diagrams have the same topology as in QCD. Moreover there are no soft but only collinear 
divergences and so factorization become simpler to deal with. For this reasons it has been 
customary, in the past, to face a QCD problem by studying it first in (0^)6 0, H, % |10|) HI 



In U-^) 6 one can define as in QCD a scalar parton density and a fragmentation function 



which obey DGLAP evolution equations ||^, 11|. It is in this framework that we are going 
to calculate here the one particle deep inelastic cross section in the limit t <^ Q^. The aim 
of such a calculation is twofold. On one hand we want to verify that the renormalized one 
particle cross section gets large logQ'^/t corrections, as one naturally expects in the two 



1 



scale regime t <^ Q^. On the other hand we want to check the consistence of our result 



with the factorization put forward in Ref.[^ 



This paper is organized as follows: in Sect. 2 we will recall the result one gets in (0^)6 
for inclusive DIS. In Sect. 3 we will perform the calculation of the semi-inclusive cross 
section. In Sect. 4 we will show how an interpretation in terms of cut vertices can be 
given. In Sect. 5 we will sketch our conclusions. 



2 DIS 



As a first step we will focus on the deep inelastic inclusive cross section. We consider the 
process p + J(g) X where J = is a scalar operator which plays the role of the 
electromagnetic current. We define as usual 



^2 2 

Q = -q x = — . 

2pq 



The structure function is 

12 



f2.11 



W{x,Q') = d%e'^y <p\J{y)Jmp> . (2.2) 

We want to calculate the parton-current cross section w{x, Q^) at = in dimensional 



Figure 1: Lowest order contribution to the deep inelastic cross section 
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regularization {D = 6 — 2e), by retaining only divergent contributions. At Born level we 
have (see Fig. 1) 

Wo{x, Q2) = ^2n6{{p + qf) = 5(1 - x). (2.3) 
The first order corrections are shown in Fig. 2 and give 

w,a{x, Q') = 7^x(l - x) (--) (^) (2.4) 



(47r)3 ^ ' \ ej \ Q 
-i^(-^Q^)-2xl^l(lp!)%(l-.) (2.5) 

wUx, Q^) = 1-^ (--] ( ^] ^(1 - X). (2.6) 



12 (47r)3 V e J \ 

The factor 2 in eq. ( |2.5| ) takes into account the contribution from the symmetric diagram. 
The diagram in Fig. 2d gives only a finite correction. External self energies have not 
be taken into account^] since we work at p"^ = 0. Usually one has not to worry about 
current renormalization, because the electromagnetic current is not renormalized by strong 
interactions. However in our simplified model the operator Jq = j4>o is renormalized by 
the interaction. The renormalization constant Zj is defined by 

J = ZjVo. (2.7) 

We get, in the MS scheme 

As a matter of fact, being the renormalization scale dependent, the coupling to the cur- 
rent becomes scale dependent, so it turns out more convenient to define a dependent 
renormalization 



12 (47r)3 e V Q 



^ Actually one should also include self energy corrections on the current line, which are of the same 
order in A and cause the mixing of the operator cf)'^ with the operator Nevertheless these diagrams 

are vacuum polarization of the external probe, and in the real world are suppressed by one power of aem 
and so are neglected. We decide here to follow the literature |l^, |l) and not to include them in our 
definition of the structure function. 
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(a) 



(b) 



(c) 



(d) 



Figure 2: One loop corrections to the deep inelastic cross section 

which takes into account this effect. Therefore the cross section is obtained summing up 
all the contributions and multiplying by Zj^(Q^). Up to finite corrections we get 



wix,Q^) =5il-x) + 



{Any 



:P(x) 



2\ <: 



(2.10) 



where 



P{x) = x{l-x) - —S{l-x) 



(2.11) 



is the DGLAP kernel for this model [|, |T4|. We see that this result has the same structure 
one gets in QCD. The contribution to the structure function is obtained as a convolution 
with the bare parton density fo{x) 

du 



W{x,Q' 



u 



-foiu)w{x/u,Q ). 



(2.12) 
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The collinear divergence in w{x, Q^) can be absorbed as usual by defining a dependent 
scalar parton density f{x, Q"^) by means of the equation 



Jx U 



(5(1 -m) + — — P(m)- ' ^ 



fix/u,Q'). (2.13) 



(47r)3 ' 'e \ 

This renormalized scalar parton density obeys the DGLAP evolution equation 

Q'^f(^^Q')= / -P{u)f{x/u,Q'). (2.14) 

On the same footing one can define for the process J{q) ^ p + X with q timelike 
a fragmentation function d{x, Q"^) which obeys the same DGLAP evolution equation. 
Thanks to the Gribov-Lipatov reciprocity relation at one loop level the timelike 
DGLAP kernel is the same as in the spacelike case, but this relation is broken at two 
loop, as explicitly verified for (0^)6 in Ref.p. 



3 Semi-inclusive DIS 

We are now ready to discuss the process p + J{q) p' + X. In this case a semi-inclusive 
structure function can be defined as 

■)2 



W{p,p',q) = |-E / ^'^e'"' <p\J{x)\p'X><Xp'\J{0)\p 



> . (3.1) 



X 

We define 

. = ^. (3.2) 
pq 

We will calculate the partonic cross section in the limit 

t < g2 (3 3) 

at leading power, by keeping only divergent terms and possible logQ^/t contributions. It 
turns out that the approximation ( p.3|) selects a special class of diagrams, those in which 
the produced particle is radiated by the incoming one. As a matter of fact at the lowest 
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Figure 3: Leading order contribution to the one particle deep inelastic cross section in the 
region t <^ 

order in A there is no contribution in the region t <^ Q^. The first diagram is the one in 
Fig. 3 which gives 



where Ao is the bare couphng constant. Notice that, at this order it is necessary to 
distinguish between A and Aq. 

The other contributions come from the one loop corrections hsted in Figs. 4,5. It 
appears that the topology of the diagrams is much richer than the one for the inclusive 
scattering. The diagram in Fig. 4(a) gives 



xS{l — X — z) 



(3.4) 



(3.5) 



We choose as a pair of lightlike vectors p and q' = 



q + xp and set 



k — ap + (3q + ]s.± 



(3.6) 



p' = a'p + (3'q + p'^ 



(3.7) 



s — 2pq. 



(3.8) 
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Figure 4: One loop leading contributions to the one particle deep inelastic cross section 
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Figure 5: One loop finite corrections to the one particle deep inelastic cross section 

The delta functions can be used to integrate over a and (5. It turns out that the integral 
over kj^ gets the leading contribution in the region of small k^^. In this limit we have 

W2a{x, Z, t, Q'^) ~ ^X{1 -X- zfX^ 

J (27r)^-i ki ((1 - z)kl - 2(1 - X - z)p'^k^ + {l-x-z){x + z)tf ^ ' ' 
so there is a collinear divergence at kj^ = which corresponds to the configuration in 
which k becomes parallel to p. The result is 

In the same way the diagram in Fig. 4(b) has a collinear pole when k gets parallel to p' 
and we get 

= j^^x^^ (--j (4^j . (3.11) 
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The other calculations are straightforward and give 

^.(.,.,t,Q^)^^^^(l-^)log^ (3.12) 

w,,{x, t, Q') = l-^^x (--] (^] 5(1 - X - z) (3.13) 



12 (47r)3 t2 \ ej \ Q 



W2eix, z, t, Q2) ^ 2 X i^^x (^^^ 5(1 - X - ;.) (3.14) 



w,A^^ Q') = 2 X ^^^x Q 5(1 - X - .) (3.15) 



W2,{x,z,t,Q^) = 2^—-——x[--\[^ \ b{\-x-z). (3.16) 



12 (47r)3 ^2 \^ e y V ^ 



Again the factor 2 in eqs. (p.l4|) - (|3.16|) takes into account the contribution of the sym- 
metric diagrams. The diagrams in Fig. 5 are finite and don't give logQ^/t contributions 
so they will be neglected. 

There are of course many other diagrams which have not been considered here. They 
are either diagrams in which the produced particle comes from the fragmentation of the 
current and diagrams which represent interference between target and current fragmenta- 
tion. However one can verify immediately that these diagrams are suppressed by powers 
of t/Q^, and so they are not relevant here. For example, consider the diagram depicted in 
Fig. 6(a), in which the observed final state particle comes from current fragmentation. The 
explicit calculation shows that this diagram is suppressed by two powers of tjQ^ . This 
fact can be better understood if we redraw it as in Fig. 6(b): in the large limit there 
are actually four lines going into a hard vertex and the diagram is suppressed by power 
counting ||12|. Summing up all the contributions, multiplying by Zj^((5^), introducing 
the running coupling constant 

= l + /3oA4ogt/4vr/.2 - - i (ip^' (3.17) 
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(a) 



(b) 



Figure 6: (a) A current fragmentation contribution to the semi-inclusive cross section 
(b) The same diagram depicted in a different form 

and using 
we finally get 



w{x,z,t,Q ) = ——xi6{l - X - z) + -6{1 - X - z) 

t^ \ 6 [AiTy e \ t 

/ 1 — X — z 1 — X — z\ / 1\ /47ru^^ " 



(47r)3 V (x + 2)2 (l-xyj\ ej\ t 



X \1- zj (47r)3 t 
We have now to verify that the collinear divergences in eq. (|3.19| ) can be correctly absorbed 



in the redefinition of parton density and fragmentation function. We have 

du dv 

W{x,z,t,Q^) = — —fo{u)w{up,p'/v,q)do{v) 
Jo U Jq V 

— f —rfo(u)w(x/u,z/uv,tu/v,Q'^)do(v). (3.20) 

u Jo 

Here the integration measure dv/v^ takes into account the scaling of the phase space 
element of the detected particle. 
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By using the definition of scale dependent parton density ( |2.13| ) and tlie corresponding 
definition for fragmentation function, eventually we find 



W{x,z,t,Q' 



^ du dv 



x+z 



- / -^f{u,t) 



u 



6 1 



x/u 



1 — z/uv 



)i2 



x/u 



(47r)^ \l — z/uvj ^ t - 



d{v, t) 



(3.21) 



where again only leading logQ^/t terms have been considered and the integration limits 
are derived using momentum conservation. 



Eq. ( |3.21| ) is our final result. We see that the renormalized hard cross section gets, 
as expected, a large logQ^/t corrections whose coefficient is precisely the DGLAP kernel. 
Such correction, if not properly resummed, can spoil perturbative calculations in the 
region t <^Q^ . 

Moreover we find that in the limit t — a new singularity appears in the cross section, 
which corresponds to the configuration in which p becomes parallel to p' . As pointed out 
in Ref. 0], when integrating over t, such singularity can't be absorbed in ordinary parton 
densities and fragmentation functions and the introduction of a new phenomenological 
distribution, the fracture function becomes necessary. 



Eq. ( ^.21| ) can also be rewritten in the following form (see Fig. 7): 

r 



W{x,z,t,Q'' 



zt^ 



+ T^Pf-)l0g 



-p(-)f{u,t) 

^^j. u[u — r) \u 



6 1- 



X 



(Att)^ \rJ "tJ \u — r 



d 



z 



t 



(3.22) 



where we have defined the A-P real scalar vertex P{x) = x{l — x). The function 
E^'^ (x, Q^ Ql) = 6{1 -x) + j^,P{x) log % 



(3.23) 



is the first order approximation of the evolution kernel E{x,Q^,Qq) which resums the 
leading logarithmic series [14 . This fact suggests that an interpretation of eq. ( |3.22 



can 



be given in terms of Jet Calculus |T^ 
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Figure 7: Graphical interpretation of eq. ^ 



4 Interpretation in terms of cut vertices 



The cut vertex expansion is a generalization of the Wilson expansion originally proposed 



by Mueller in Ref .11161 and applied to a variety of hard processes in Ref.|p.7||. We will 



briefly recall it for DIS in (0^)6- 

Let us go back to Sect. 2 and set < with p = 0,p^). If we choose a frame in 



which ^ p_ we can write for the parton-current cross section |]T6 
du 



w{p, q) 



u 



-v{p\u)C{x/u,Q') 



where 



j V{p, k)x5 — 



(4.1) 



(4.2) 



P+J (2vr)6 

is a spacelike cut vertex, which fully contains the mass singularities of the cross section, 
and C(x,Q^) is the corresponding coefficient function. In eq. (|4.2|) V{p,k) is defined as 
the discontinuity of the four point amplitude in the channel {p — kY and the integration 
over k is properly renormalized. 
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The results of the previous sections can be easily recast in the form of cut vertex 
expansion. If we define 

v{x,e) = 5{l-x) + -^Pix)(--] (4.3) 



{Any 

Q') = 6(1 -x) + ^P(x) log ^ (4.4) 
we can write eq. ( |2.10| ) in the form 

w(a;,Q2)=/ —v{u,e)C{x/u,Q^). (4.5) 

Jx U 

Here v{x, e) is a spacelike cut vertex defined at = whose mass divergence is regularized 
dimensionally. 




Figure 8: Factorization of the semi-inclusive structure function 
It turns out that a similar interpretation can be given of eq. ( |3.19| ). We define 



X 



X 



1 - z 



(4.6) 



and 



v{x, z, t, e) 



t2 



6{l-x) + P{x)j^\og—^+ ' 



t ■ (47r)3 e V t 
1 — z)^x(l — X 



X i'sil X) (l-^)^^(l-^) 

V6^ ^ (a;(l-z) + z)2 [i - x{l - z)y 



(4.7) 



13 



as a generalized cut vertex |]T2| which contains all the leading mass singularities of the 
cross section. We can write up to 0{t/Q'^) corrections (see Fig. 8) 

w{x,z,t,Q'^,e) = / —v{u,z,t,e)C{x/u,Q^) (4.8) 

where the coefficient function is the same which occurs in inclusive DIS. 

It is well known that in dimensional regularization there is a mixing between collinear 
and ultraviolet divergences. In order to avoid it one should distinguish between e and 
e' to regulate UV and collinear divergences, respectively. Moreover external self energies 
should be taken into account since they are zero on shell for a cancelation of the two kind 
of divergences. In this framework one can show that factorization ( [4.8| ) actually holds 
graph by graph [[T^ . 



5 Conclusions 

In this paper we have studied the deep inelastic semi-inclusive cross section in the target 
fragmentation region and we have performed an explicit calculation in (0^)6 model field 
theory. We have shown that the renormalized hard cross section gets a large log Q'^/t 
correction as expected in a two scale regime. Furthermore we have found that the coeffi- 
cient driving this logarithmic correction is precisely the scalar DGLAP kernel. This result 
suggests that the dependence of the cross section in such processes at fixed z and t is 
driven by the same anomalous dimension which controls the inclusive DIS, as proposed 
in Ref. [|r^, and in in the context of diffraction in Ref.[0. 

We have then examined our result from the point of view of extended factorization 
and we have found that it is consistent with such an hypothesis. In this framework 
the partonic semi-inclusive cross section factorizes into a convolution of a new object. 



a generalized cut vertex v{p,p',x) ||T2|, with four rather than two external legs, and a 
coefficient function C{x,Q^). The former is of long distance nature and embodies the 
leading mass singularities of the cross section, while the latter is of perturbative nature 
and, what is important, it is the same as in inclusive DIS. 
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Therefore these results verify the vahdity of the approach proposed in Ref. IT^. Of 
course the calculation performed here is only a one loop calculation in a scalar model. 
Nevertheless we believe that the results obtained maintain the same structure in QCD, 
by assuming in (0^)6 a particularly simple and appealing form. 
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